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CHAPTER  I 


INTRODUCTION 

1. 1  General  Introduction 

Elastic  contact  problems  and  elastic  indentation  problems 
involving  the  normal  application  of  a  load  against  an  elastic  half¬ 
space  by  a  rigid  body  have  been  of  considerable  interest  in  various 
fields  of  applied  mechanics.  The  point  indentation  problem  was 
solved  in  closed-form  by  Boussinesq  [1]  and  Michell  [2].  The  closed- 
form  solution  for  spherical  indentation  was  developed  by  Hertz  [3] 
and  Huber  [4].  The  conical  indentation  problem  was  solved  in  closed- 
form  by  Sneddon  [5,6]. 

As  the  theory  of  plasticity  has  advanced,  elastic/plastic 
contact-indentation  problems  have  been  approached  analytically  [7,8,9], 
experimentally  [10,11]  and  numerically  [12,13,14,15].  To  the  author's 
knowledge,  the  problem  of  normal  loading  of  an  elastic/perfectly 
plastic  half-space  by  a  rigid  conical  indenter  has  not  been  analyzed, 
despite  the  fact  that  aspects  of  this  problem  are  of  great  interest 
in  many  fields,  including  indentation  fracture  mechanics.  An  analytical 
solution  to  this  problem  would  be  extremely  difficult  to  develop  but 
the  use  of  the  high-speed  computer  combined  with  the  proper  finite 
element  code  have  made  its  solution  possible.  In  this  investigation  a 
finite  element  code  which  incorporates  the  pure  tangential  stiffness 
incremental  method  was  applied  to  obtain  a  solution  to  the  problem  of 
elastic/plastic  conical  indentation. 


The  purpose  of  this  investigation  was  to  perform  fundamental 
research  on  elastic/plastic  conical  indentation  problems  associated 
with  the  machining  of  brittle  materials.  It  is  anticipated  that  the 
design,  fabrication  .and  maintenance  capabilities  for  brittle  structures 
can  be  improved  when  the  basic  understanding  of  fundamental  failure 
mechanisms  associated  with  resultant  generated  elastic/plastic  stress 
fields  and  residual  stress  fields  is  increased.  The  resultant 
axisymmetric  elastic/plastic  stress  fields  associated  with  both 
loading  and  unloading  can  then  be  related  to  median,  radial  and 
lateral  cracking.  In  order  to  obtain  this  solution,  it  was  necessary 
to  apply  axisymmetric  elastic/plastic  analysis  and  cyclic  loading 
in  the  finite  element  program. 


1.3  Scope  of  the  Investigation 

This  investigation  consists  of  four  primary  sections.  The  first 
section  considers  the  development  of  a  finite  element  model  for  the 
conical  indentation  problem  and  the  verification  of  the  model  by 
comparing  numerically  generated  results  with  the  closed-form  results 
for  an  elastic  half-space  derived  by  the  well-known  Sneddon  solution. 
Secondly,  an  elastic/plastic  conical  indentation  analysis  is  performed 
using  the  finite  element  model  with  plastic  analysis  capability  to 
generate  the  elastic/plastic  stress  and  strain  fields  and  to  determine 
the  extent  of  plastic  deformation  near  the  indentation  site.  Thirdly, 
a  loading-unloading  cycle  is  simulated  in  order  to  obtain  the  residual 


stress  fields  associated  with  elastic/plastic  indentation.  Finally, 
the  results  in  terms  of  generated  stress  fields  are  related  to  the 
initiation  and  propagation  of  median,  radial  and  lateral  cracks 
beneath  the  indenter. 


CHAPTER  II 


ELASTIC/PLASTIC  FINITE  ELEMENT  METHOD 

2.1  General  Concepts 

This  section  discusses  the  general  concepts  utilized  in  the 
solution  of  elastic/plastic  problems  hy  finite  element  techniques. 

This  approach  is  based  on  a  plastic  stress-strain  matrix  derived 
from  the  Prandtl-Reuss  equations,  and  on  incremental  or  iterative 
stiffness  methods. 

2.1.1  Stress-Strain  Equations 
It  has  been  well  established  that  the  plastic  constitutive 

equation  can  be  expressed  in  incremental  form.  In  this  investigation, 
the  selected  incremental  plasticity  relations  conform  to  the  Prandtl- 
Reuss  flow  rule,  which  is  the  usual  flow  rule  associated  with  the 
well-known  Huber-Mises  yield  surface  criterion  [7,16,17].  The 

Prandtl-Reuss  flow  rule  uses  the  strain  increment  de.  which  is 

i.r 

related  to  the  stress  increment  do, .  .  The  incremental  stress- 

ii 

strain  relations  with  the  differential  form  of  the  Huber-Mises 
yield  surface  criterion  can  be  represented  in  matrix  form  as: 

{do}  =  [Dp]{de>  ,  (1) 

where  {do}  and  {dt}  are  defined  as  the  .olumn  matrices  of 
do^j  and  dt_^,  respectively;  and  [DP]  is  the  plastic  stress-strain 
matrix.  Equation  (1)  can  be  expressed  as: 


(i  =  1,2, 3, 4, 5, 6), 


(2) 


6  3o . 

da  =  I  T-i  de. 
i  .  ,  3e.  j 
J=1  J 


where  the  coefficients  30^3^  are  the  partial  stiffness  components 
[28].  Hooke's  law  for  isotropic  elastic  materials  can  be  expressed  in 
matrix  form  as: 


{a}  =  [De] { e }  , 


(3) 


where  the  elastic  stress-strain  matrix  [D  ]  is  given  as: 
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For  the  elastic/plastic  finite  element  method,  the  plastic  stress- 
strain  matrix  [D^]  for  yielded  elements  takes  the  place  of  the  elastic 
stress-strain  matrix  [D  ]. 

The  Prandtl-Reuss  stress-strain  relations  for  the  deviatoric 


strain  increment  de^.  during  continued  loading  are: 
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The  Huber-Mises  yield  criterion  and  its  differential  form  are: 
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From  equations  (5)  and  (11) , 
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The  total  stress  increment  do..  is: 

ij 


da . .  =  da !  3  + 


6  . de  .  . 


ij  uij  3(1- 2v)  "ij  ii 


m 


V 


=  2G 


.  .  v  „  ,  ,  °k£dekl 

[deij  +I^  5ijdeii  •  aij  — S 


Equation  (15)  can  be  represented  in  matrix  form  as: 


where 
[Dp]  «  2G 
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It  is  obvious  that  equation  (16)  is  similar  in  form  to  equation  (3)  and 


that  the  necessary  modifcation  from  elastic  to  elastic/plastic  behavior 
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0  n 

would  be  to  replace  [D  ]  by  [D  ]  for  yielded  elements  as  successive 
yielding  takes  place.  This  procedure  was  developed  by  Marc3l  and 
King  [18]  and  Yamada,  Yoshimara  and  Sakurai  [19,20]. 

2.1.2  Incremental  Stiffness  Method 

In  this  investigation,  the  incremental  pure  tangential  stiffness 
method  was  applied  to  construct  the  stiffness  matrix.  In  this  section, 
this  incremental  stiffness  method  is  discussed.  For  elastic  finite 
element  analyses,  the  relationship  between  the  nodal  force  and  nodal 
displacement  in  each  element  n  can  be  expressed  in  matrix  form  by: 

(Q>n  -  [k6]  {q>n  ,  (18) 

where 

[ke]  -  ///  [B]T[De][B]dV  ,  (19) 

and 

<G>n  *  [B]  {q>n  ,  (20) 

(Q)n  =  equivalent  nodal  force  matrix, 

(q>n  “  equivalent  nodal  displacement  matrix, 
[B]  =  displacement-strain  matrix. 

Integration  of  equation  (19)  is  over  the  volume  of  the  element.  The 
matrix  [ke]  is  called  the  stiffness  matrix  for  the  elastic  elements. 
Similarly,  for  the  plastic  elements  the  stiffness  matrix  is  expressed 
as: 

[kP]  =■  fff  [B]T  [DP]  [B]  dV  , 


(21) 
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and 

{dQ}  =  [ kP] { dq }  .  (22) 

n  n 

The  overall  stiffness  matrix  [k] ,  which  is  an  assemblage  of  the  stiffness 
matrices  fk  ]  and  [k*j,  relates  the  nodal  load  increment  { dQ }  to  the 
nodal  displacement  increment  {dq}  as: 

{dQ}  =  [k]  {dq}  .  (23) 

In  order  to  solve  the  incremental  equation  (22)  for  an  elastic- 
plastic  problem,  the  load-deflection  relation  is  first  expressed  by: 

{Q}n=  [k]{q}n 

-  /v  [B]T  [D]  [B]  dV  {q}n 

=  /v  [B]T  {o }dV  .  (24) 

Using  variational  techniques, 

{dQ}n  =  /v  [B]T  {dcr}  dV 

-  /v  [B]T  [Dt1  {dt }dV 

-  /v  [B]T  [Dt3  [B]  dV  {dq} 

-  [1^]  {dq}n  ,  (25) 

where  [D^]  and  [kj]  are  called  the  material  tangential  modulus  matrix, 
and  the  structural  tangent  stiffness  matrix,  respectively.  It  is 
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obvious  that 

[D6]  »  [Dx]  , 
[ke]  =  [1^] 

for  the  elastic  element  and 


[DP]  =  [Dx]  , 

[kP]  =  [kT] 

for  the  plastic  element. 

The  matrix  [DP]  is  updated  for  each  increment  of  load  with 
computation  of  the  plastic  strain  increment  {deP}  at  the  yielded 
elements.  The  tangent  stiffness  matrix  [k^]  is  computed  at  the  end 
of  each  increment  and  used  for  each  succeeding  increment  according 
to  the  following  equations: 


(26) 

(27) 

(28) 


after  the  ic^  increment  where  (Q  }  and  (q  }  are  the  initial  loads 

o  o 

and  displacements,  usually  null  vectors.  The  above  procedure  is 
called  the  "incremental  tangent  stiffness  method." 
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2.1.3  Other  Stiffness  Methods 

For  elastic/plastic  finite  element  analyses,  there  are  three 
major  stiffness  methods  in  addition  to  the  incremental  stiffness 
method,  the  iterative  method,  the  initial  strain  method  and  the  initial 
stress  method.  The  iterative  method  [21,22]  is  used  for  elastic/plastic 
behavior  when  the  deformation  theory  of  plasticity  is  employed.  Suppose 
the  nonlinear  nodal  force  and  nodal  displacement  relation  is  expressed 
as : 

{Q}n  *  [k(q,Q)]fq}n  ,  (29) 

where  the  force  {Q}^  is  entirely  known  and  the  displacement  {q}^  is 
entirely  unknown.  Such  an  iterative  procedure  involves  an  initial  guess 
for  the  displacement  .  An  increasingly  accurate  series 

vectors 

{q}n1}’  {q}n2>’  •  *  •  {q}nN)  (30) 

is  then  generated  with  the  objective  being  convergence  to  the  exact 
vector  (q }  .  The  first  iteration  is  obtained  from: 

[k((q}^0)  ,  (Q}n)]  {q}<1}  =  {Q}n  ,  (31) 

which  has  a  formal  solution: 

{q}^  =  [k(°)]'1  {Q)n  .  (32) 

The  general  procedure  is  to  solve  the  equation: 

[k({q}(i_1),  { Q }  ) ]  {q}(i) 


=  {Q)  . 


(33) 
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In  order  to  find  the  current  step,  the  displacements  from  the  pre¬ 
ceding  step  are  used. 

The  major  problems  arising  from  the  use  of  this  method  are: 

1)  there  is  no  guarantee  of  convergence, 

2)  a  new  secant  stiffness  must  be  generated  at  each  step,  and 

3)  a  new  stiffness  matrix  must  be  inverted  at  each  step. 

The  preceding  two  methods,  incremental  method  and  iterative  method, 
are  variable  stiffness  methods  with  stiffness  matrices  generated  for 
each  step.  In  addition  to  these  variable  stiffness  methods,  there  are 
constant  stiffness  methods:  the  initial  strain  method  and  the  initial 
stress  method.  These  have  been  developed  to  reduce  the  computation 
time  by  utilizing  the  same  stiffness  throughout. 

In  the  case  of  the  initial  strain  method  [22],  use  is  made  in  the 
first  iteration  of  the  elastic  stress-strain  [De]  given  by  equation  (4) 
to  define  [k  ]  of  equation  (19)  and  solve  the  equation: 

[ke]  (q}(1)  =  {Q}  ,  (34) 

where  {q}^  defines  the  displacement  vector  for  the  first  iteration 
and  (Q }  is  the  total  prescribed  load  vector.  This  yields  an  approxi¬ 
mate  displacement  vector  {q}^.  A  corresponding  approximate  strain 
is  found  from: 

(e}(1'  =  [B]  {q }  (1> .  (35) 

Hence  the  exact  stress  corresponding  to  the  approximate  strain  is: 

{o}(1)  =  [De]  |{c}(1)  -(eo}(1)}  , 


(36) 
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where  {e  is  an  artificially  imnosed  prestrain.  An  artificial 

o  ' 

equivalent  nodal  load  corresponding  to  (e  is: 


{Q}^  ,  -  A,  [B]T  [De]  {e  }(1)  dV  . 
e , n,  l  V  o 


(37) 


The  second  iteration  becomes: 


[ke]  {q}(2)  =  {Q}  +  {Q}^  #  , 


(38) 


(2) 

in  which  the  second  displacement  vector  {q}  is  sought.  The 
recurrent  equation  then  becomes: 


[ke]  (q } (i)  =  (Q)  +  . 


(39) 


For  the  initial  stress  method  [22,23,24],  equation  (36)  can 
be  described  as: 


{a} (1)  =  [De]  {t}(1)  +  (oo>  ,  (40) 

where  {a  is  an  artificially  introduced  prestress.  Corresponding 

to  {a  }  we  have 
o 


(1)  =  . 
e,n,2. 
o„ 


A,  [B] 


{<V 


(1) 


dV 


The  second  iteration  involves  solution  {q} 


(2) 


[k6]  (q}(2)  =  {QJ  +  . 

e ,  n ,  4 


a 

o 


as : 


(41) 


(42) 


The  recurrent  equation  is 


[k6]  {q}(i)  =  {Q}  +  (43 

e,n ,  £ 

o 

o 


2.2  BOPACE  -  3D 

Selection  of  the  proper  finite  element  program  that  has  plastic 
analysis  capability  is  an  important  aspect  of  solving  the  elastic/ 
plastic  indentation  problem.  Here,  a  finite  element  code  BOPACE  (THE 
BOEING  PLASTIC  ANALYSIS  CAPABILITY  FOR  ENGINES)  -3D  Version  6.2,  which 
was  developed  for  analysis  of  the  Space- Shuttle  Main  Engine  (SSME) , 
was  chosen. 

The  program  employs  the  Huber-Mises  yield  surface  criterion,  the 
Prandtl-Reuss  flow  rule,  and  an  isotropic-kinematic  hardening  theory 
[7,25,26,27]  as  the  basic  concepts  of  the  elastic/plastic  theory.  It 
also  employs  three  types  of  stiffness  methods  [28], 

1)  the  pure  "tangent  stiffness  incremental"  method, 

2)  the  "constant-stiffness,  initial  stress"  method,  and 

3)  the  "combined"  method 

for  solution  of  elastic/plastic  problems,  and  a  modified  Gauss  wave- 
front  solution  procedure  [22]  for  linear  equation  solution. 

Development  of  the  program  has  been  strongly  influenced  by  the 
following  requirements: 

1)  the  analysis  of  very  high  temperature,  large  plastic- 
creep  effects  and  geometric  nonlinearities, 

2)  treatment  of  cyclic  thermal  and  mechanical  loads. 


3)  improved  material  constitutive  theory  which  closely  follows 
actual  behavior  under  variable  temperature  conditions,  and 

4)  a  stable  numerical  solution  approach  which  avoids  cumulative 
errors . 

The  above  characteristics  make  it  one  of  the  most  sophisticated  finite 
element  codes  currently  available  for  many  types  of  general  nonlinear 
problems. 

2.3  Verification  of  Nonlinear  Finite  Element  Code 

Before  applying  the  solution  technique  adopted  in  this  investiga¬ 
tion  to  the  indentation  problem,  the  technique  was  first  verified 
by  comparing  generated  elastic/plastic  results  for  two  variations  of 
a  selected  plane-stress  problem  with  those  independently  obtained  from 
the  literature.  The  plane-stress  problem  selected  was  a  plate  with  a 
circular  penetration  subjected  to  uniaxial  tension. 

2.3.1  Comparison  to  Yamada' s  Analysis 

Here,  the  results  obtained  from  the  elastic/plastic  finite  element 
code  used  in  this  investigation  are  compared  to  those  obtained  by 
Yamada  [Personal  Note]  for  an  aluminum  plate  with  a  circular  penetration 
subjected  to  uniaxial  tension.  The  comparison  is  made  in  terms  of 
forces  and  displacements  along  the  vertical  and  horizontal  symmetry 
axes  of  the  plate  for  the  same  nodal  coordinates.  The  Yamada  solution 
employs  an  incremental  theory  upon  which  many  current  nonlinear  finite 
element  codes  are  based  [NASTRAN,  NONSAP,  etc.] 


The  stress-strain  curve  and  pertinent  material  properties  for 
the  aluminum  plate  material  are  shown  in  Figure  1  for  both  isotropic 
and  kinematic  hardening.  Figures  2  and  3  show  the  finite  element 
meshes  adopted  by  Yamada  and  in  this  investigation,  respectively .  Both 
meshes  employ  the  same  nodal  point  coordinates  for  easy  comparison  of 
solutions  but  different  element  shapes  since  Yamada' s  program  considers 
only  triangular  elements  and  the  finite  element  code  adopted  in  this 
investigation  utilizes  quadrilateral  elements.  This  resulted  in  an 
apparent  distortion  of  quadrilateral  elements  (Figure  3).  Due  to 
symmetry,  only  the  first  quadrant  of  the  plate  was  discretized.  A 
uniform  tensile  stress  of  39.2  N/mm^,  which  was  selected  in  the 
Yamada  solution,  was  applied  in  the  y-direction  along  the  upper  edge 
of  the  plate. 

Comparisons  of  the  nodal  forces  and  displacements  along  the 
horizontal  and  vertical  axes  of  symmetry  are  shown  in  Figures  4  through 
7.  Results  obtained  from  the  two  finite  element  codes  are  within  10 
percent,  with  the  differences  caused  primarily  by  the  dissimilarity 
of  element  shapes.  Theoretically,  the  accuracy  could  be  improved  for 
a  given  number  of  nodes  by  using  quadrilateral  elements,  since  the 
increased  degrees-of-freedom  would  permit  a  closer  approximation  to 
the  displacements  within  an  element  [22,23]. 

2.3.2  Comparison  to  Marcal  and  King's  Analysis 


In  the  previous  section,  a  coarse  grid  model  was  used  to  compare 
two  different  solution  techniques.  In  order  to  determine  the  extent 
and  character  of  the  plastic  region  as  a  function  of  applied  tensile 
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Figure  1.  Stress-Strain  Curve  and  Pertinent  Material 
Properties  for  Yanada's  Analysis. 
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Figure  2 


Yamada's  Finite  Element  Model  for  Aluminum 
Plate  with  Circular  Penetration  under 
Uniaxial  Tension. 
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Figure  3.  Finite  Element  Model  Used  in  Thi 
Investigation . 
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Figure  5.  Nodal  Forces  F  along  the  y-axis. 
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Figure  6.  Nodal  Displacement  U,  along  the 
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Figure  7.  Modal  Displacements  U  along  the  y-axis. 
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stress,  the  model  was  reconstructed  so  that  nodal  density  near 
the  stress  concentration  would  be  increased.  Through  this  remesh¬ 
ing,  the  number  of  nodes  and  elements  and  the  maximum  wavefronts  were 
increased  to  399,  350  and  33  from  94,  70  and  15,  respectively.  This 
improved  model  is  shown  in  Figure  8.  The  material  used  for  this 

investigation  was  linear  strain-hardening  aluminum,  with  the  following 

2 

mechanical  properties;  Young's  modulus  7000  kg/mm  ,  Poisson's  ratio 

2  2 
0.2,  Yield  stress  24.3  kg/mm  ,  strain  hardening  rate  225  kg/mm  . 

These  values  were  obtained  from  Zienkievicz 's  [21]  and  Marcal  and  King's 

[18]  finite  element  analysis  of  a  linear  strain-hardening  aluminum 

perforated  strip,  whose  diameter-width-height  ratio  was  1:2: 3. 6 

instead  of  1:3:5  of  this  investigation. 

The  plastic  regions  for  load  factors  of  1.5,  2.0  and  2.5  are 

shown  in  Figure  9.  A  load  factor  of  1.0  represents  the  load  required 

to  cause  incipience  of  yielding  in  the  plate.  Load  factors  of  1.5, 

2.0  and  2.5  resulted  in  progressive  yielding  of  the  plate  near  the 

circular  penetration.  Figure  10  shows  the  plastic  regions  obtained 

by  Marcal  and  King.  Comparison  of  Figure  9  and  10  indicates  that  the 

character  and  extent  of  the  plastic  regions  are  nearly  identical 

for  both  analyses  in  spite  of  the  shape  difference  of  the  specimens. 
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Figure  3.  Improved  Model  used  in  the  Investigation 
of  Plastic  Regions. 
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Figure  9.  Character  and  Extent  of  Plastic  Regions 
Obtained  in  this  Investigation  for 
Selected  Load  Factors. 


CHAPTER  III 


SOLUTION  TECHNIQUES  AND  THE  INDENTATION  MODEL 
3.1  Method  of  Approach 

Since  an  analytical  solution  to  the  elastic/plastic  problem 
is  not  available  at  this  time,  the  investigation  of  elastic/plastic 
conical  indentation  was  accomplished  entirely  by  finite  element 
techniques.  The  objective  was  to  develop  and  demonstrate  numerical 
techniques  which  can  be  applied  to  elastic/plastic  problems 
involving  flow,  fracture  and  residual  stress  effects  beneath  a 
conical  indenter.  These  techniques  will  be  explained  in  greater 
detail  in  subsequent  chapters. 

In  the  loading  procedure,  three  separate  stages  were  required 
to  simulate  elastic/plastic  loading  and  unloading  of  a  half-space 
by  a  rigid  conical  indenter.  In  the  first  stage,  initial  loading 
was  applied  through  nodal  displacements  selected  to  simulate 
conical  indentation  at  the  contact-surface.  Displacement  loading 
was  used  since  there  is  no  analytical  method  currently  available 
for  predicting  the  load  distribution  for  elastic/plastic  conical 
indentation.  The  second  stage  was  designed  to  obtain  the  equivalent 
surface  nodal  loads  which  would  result  in  the  same  stress  and  displace¬ 
ment  fields  as  those  obtained  by  the  displacement  loading  of  the  first 
stage.  This  was  accomplished  by  using  the  numerical  results  obtained 
following  the  displacement  loading  of  the  first  stage.  In  the  third 
stage,  the  equivalent  surface  nodal  loads  were  reduced  to  zero  in  order 
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to  simulate  unloading  of  the  specimen.  During  this  stage,  the  re¬ 
sidual  displacement  beneath  the  indenter  and  the  residual  stress 
fields  were  generated. 

The  error  norm  for  each  stage  in  the  analysis  was  limited  to 
within  5  percent. 

3.2  Finite  Element  Model  and  Grids 

The  axisymmetric  finite  element  idealization  for  the  indentation 
analysis  of  the  half-space  is  shown  in  Figures  11  through  14.  Figure 
11  shows  the  overall  outline  of  the  model  and  Figures  12,  13,  and  14 
are  sub-grids  which  must  be  inserted  into  the  appropriate  central 
sectors  of  Figures  11,  12,  and  13,  respectively,  to  construct  the 
entire  model.  There  are  4 70  nodes  and  427  elements,  and  the  maximum 
wavefront  that  governs  CPU- time  [22]  is  15.  Loading  is  applied  in 
the  vertical  direction  along  the  y-axis,  as  shown  in  the  figures.  The 
finite  element  mesh  was  designed  to  meet  the  following  objectives: 

1)  to  reduce  the  maximum  wavefront  in  order  to  minimize 
CPU- time  and  core  storage  requirements, 

2)  to  reduce  the  occurrence  of  singularities  for  elements 
in  regions  of  primary  interest  by  considering  the  shapes 
of  elements, 

3)  to  reduce  the  number  of  elements  in  regions  that  are  of 
less  interest  by  using  curved  boundary  elements  which 
include  intermediate  nodes  at  the  edge  of  the  elements, 

4)  to  limit  the  maximum  output  to  30,000  printed  lines,  which 
is  the  maximum  number  of  printed  lines  for  JOB-category  5 
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Figure  13.  Axi symmetric  Finite  Element  Sub-Grids 
for  Indentation  Analysis  (Insert  B  of 
Figure  12) . 
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Figure  14. 


Axi symmetric  Finite  Element  Sub-Grids 
for  Indentation  Analysis  (Insert  C  of 
Figure  13). 
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of  PSU-IBM  370,  by  controlling  the  total  number  of 
elements , 

5)  to  generate  or  preprocess  the  input  data  easily. 

3. 3  Material  Used  in  this  Investigation 

The  material  selected  for  this  analysis  was  idealized  soda- 
lime  glass.  Soda-lime  glass  is  a  brittle  material  used  in  many 
fracture  mechanics  experiments  and  its  behavior  is  well  established 
for  Vickers  pyramidal  indentation  and  spherical  indentation.  It  is 
known  that  soda-lime  glass  exhibits  elastic/perfectly  plastic  behavior 
when  subjected  to  indentation  loading.  It  was  thus  an  appropriate 
material  for  the  nonlinear  finite  elem  nt  model.  The  stress-strain 
curve  for  this  elastic/perfectly  plastic  soda-lime  glass  is  shown 
in  Figure  15  together  with  its  pertinent  material  properties. 


Young's  Modulus  =*49.0  (GPa) 
Poisson’s  Ratio  -  0.22 
Yield  Strength  *  3.92  (GPa) 


Figure  15.  Stress-Strain  Curve  and  Pertinent  Material 
Properties  of  Idealized  Soda-Lime  Glass  fo 
Indentation  Analysis. 


CHAPTER  IV 


ELASTIC  AND  ELASTIC/PLASTIC  CONICAL  INDENTATION 

4.1  Introduction 

The  finite  element  model  for  conical  indentation  of  a  half¬ 
space  was  formulated  so  that  it  could  be  applied  to  both  elastic  and 
elastic/plastic  indentations.  This  required  the  method  of  loading  to 
be  carefully  considered.  For  the  case  of  elastic  indentation,  the 
pressure  distribution,  resultant  forces  and  displacements  directly 
beneath  the  indenter  can  be  obtained  from  a  closed-form  solution  by 
Sneddon,  as  shown  in  Appendix  A.  Each  of  these  parameters  can  be 
used  as  load  input  to  the  finite  element  program  at  appropriate 
nodes  beneath  the  indenter  to  model  elastic  conical  indentation. 

As  the  material  begins  to  yield,  the  pressure  distribution,  and 
resultant  forces  beneath  the  indenter  are  completely  unknown,  while 
the  displacement  field  is  partially  specified.  In  the  case  of  the 
displacement  field,  nodal  points  on  the  contact  surface  must  conform 
to  the  profile  of  the  indenter.  For  regions  adjacent  to  the  contact 
surface,  the  displacement  of  surface  nodes  would  be  automatically 
generated  as  output  of  the  finite  element  code. 

In  order  to  generate  a  model  which  can  apply  to  both  clastic 
and  elastic/plastic  conical  indentations,  it  was  decided  to  input 
load  through  the  displacements  of  surface  nodes  along  the  contact 
surface.  This  would  allow  direct  comparison  of  the  finite  element 
results  with  those  obtained  from  the  closed-form  solution  for 
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elastic  indentation,  and  also  afford  partial  specification  of  the 
displacement  field  beneath  the  indenter  for  elastic/plastic  indenta¬ 
tion. 


4.2  Method  of  Analysis 

The  coordinate  system  adopted  in  the  analysis,  the  geometrical 
characteristics  associated  with  the  indenter  geometry,  and  the  dis¬ 
placement  loading  are  shown  in  Figure  16.  In  order  to  simulate 
Vickers  indentation,  the  included  angle  of  the  indenter  was  set  to 
a  =  136  degrees. 

An  elastic  finite  element  analysis  was  performed  initially  to 

verify  the  finite  element  model.  Here,  the  contact  surface  radius 

was  selected  as  b  *  0.30  mm,  resulting  in  a  depth  of  indentation 

v  =  (d+e)  =»  0.19  mm  that  was  computed  from  the  Sneddon  solution, 
max  r 

The  Sneddon  solution  also  allowed  specification  of  both  vertical  and 
horizontal  displacements  of  surface  nodes  within  and  adjacent  to 
the  contact  surface.  In  the  finite  element  analysis,  however,  only 
the  vertical  displacement  of  the  nodes  on  the  contact  surface  was 
specified.  Horizontal  displacement  was  computed  as  output  in  order 
to  allow  a  greater  degree  of  freedom  on  and  adjacent  to  the  contact 
surface. 

Following  the  elastic  finite  element  analysis,  the  elastic/plastic 
finite  element  analysis  w~j  conducted.  As  described  earlier,  the  dis¬ 
placement  field  beneath  the  indenter  is  only  partially  described  in 
this  case.  A  trial  elastic/plastic  run  was  made,  in  which  the 
initial  vertical  displacements  of  nodes  within  the  contact  surface 


were  obtained  from  the  Sneddon  solution  for  assumed  elastic/perfectly 
plastic  behavior  of  soda-lime  glass.  Following  this  step,  nodal 
force  resultants  were  checked  for  the  nodes  adjacent  to  the  contact 
surface  to  determine  whether  any  vertical  tensile  tractions  existed. 
This  condition,  when  encountered,  indicated  a  divergence  of  elastic/ 
plastic  surface  deformations  from  the  elastic  Sneddon  displacements, 
resulting  in  "overlap"  of  vertical  surface  displacements  with  the 
indenter  profile.  This  was  alleviated  through  trial  reductions  of 
the  indentation  depth  that  were  repeated  until  the  vertical  tensile 
tractions  were  reduced  to  zero  or  assumed  small  compressive  values. 
For  soda-lime  glass,  it  was  found  that  closed-form  computation  of 
vertical  displacements  on  the  contact  surface  resulted  in  no 
tensile  tractions  for  the  selected  contact  surface  radius  and  there¬ 
fore  applied  for  elastic/plastic  indentation. 

4. 3  Results  and  Discussions 

The  elastic  and  elastic/plastic  results  expressed  in  terms  of 
surface  stresses  and  stresses  directly  beneath  the  conical  indenter 
are  shown  in  Figures  17  to  20.  In  these  figures,  elastic  results 
obtained  from  the  closed-form  Sneddon  solution  and  the  elastic 
finite  element  analysis  are  compared  to  those  obtained  in  the 
elastic/plastic  finite  element  analysis.  In  each  case  stresses 
are  non-dimens ionalized  with  respect  to  the  factor  EtanS  =  (eE/b), 
while  distances  along  the  free  surface  of  the  half-space  and  beneath 
the  indenter  are  nondimensionalized  with  respect  to  the  contact  sur¬ 


face  radius,  b. 


Sneddon  Analysis 

Elastic  Finite  Element  Analysis 

Elastic/Plastic  Finite  Element  Analysis 


Figure  17. 


Non-dimensionalized  Surface  Stress  a 
along  the  x-axis.  xx 
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o  =  Sneddon  Analysis 

x  =  Elastic  Finite  Element  Analysis 

A  *  Elastic/Plastic  Finite  Element  Analysis 


Figure  19.  Non-dimensionalized  Subsurface  Stress  s  / EtanB 
and  a _ /EtanS  along  the  y-axis. 
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4.3.1  Elastic  Conical  Indentation 


Considering  initially  the  elastic  conical  indentation,  it  is 
obvious  that  excellent  agreement  is  obtained  between  elastic  stresses 
generated  by  the  closed-form  Sneddon  solution  and  the  elastic  finite 
element  analysis  both  on  the  surface  of  the  half-space  and  directly 
beneath  the  conical  indenter.  Minor  differences  are  found  at  the 
edge  of  contact  (x/b  =  1.00).  These  differences  are  caused  by  the 
inability  of  the  finite  element  grid  to  actually  express  the  high 
stress  gradient  in  this  region.  Alon<*  the  surface  at  the  edge  of 
contact,  both  the  Sneddon  solution  and  the  elastic  finite  element 
analysis  show  a  tensile  peak  in  what  will  be  termed  the  horizontal 
stress  component,  o  ,  Figure  17.  This  peak,  also  found  in  spherical 
indentation,  is  usually  associated  with  ring  cracking  and  subsequent 
Hertzian  cone  cracking.  This  will  be  discussed  further  in  a  sub¬ 
sequent  section.  Within  the  contact  zone,  this  stress  becomes 
compressive,  reaching  a  compressive  maximum  directly  beneath  the 
indenter.  Outside  the  contact  area,  this  stress  decreases  from  the 
tensile  peak  to  zero  as  distance  increases  from  the  contact  surface. 
The  out-of-plane  hoop  stress,  ozzt  along  the  surface  is  everywhere 
compressive,  with  the  peak  compressive  value  being  found  directly 
beneath  the  indenter.  Figure  18.  This  stress  also  gradually  decreases 
to  zero  as  distance  from  the  origin  increases. 

Directly  beneath  the  indenter,  the  horizontal  stress,  axx. 

Figure  19,  changes  from  high  compressive  values  at  the  surface  to 
small  tensile  values  at  approximately  one  contact  radius  from  the 


surface.  The  highest  tensile  value  is  found  at  the  point  of  y/b  =  1.60. 
Beyond  this  point,  this  stress  gradually  decreases  to  zero.  Tensile 
values  of  this  stress  result  in  median  cracking,  found  beneath 
pyramidal  and  conical  indenters.  The  subject  of  median  cracking  will 
be  discussed  in  a  subsequent  section.  The  out-of-plane  hoop  stress, 
a  ,  is  of  course  equal  to  the  horizontal  stress,  a  ,  along  this  axis 
of  symmetry.  The  subsurface  vertical  stress  component,  cr  ,  shown  in 
Figure  20,  is  everywhere  compressive;  reaching  a  peak  magnitude  directly 
beneath  the  indenter  and  rapidly  decreasing  to  zero  at  increased 
distance  from  the  origin. 

Figure  21  shows  the  pressure  distribution,  o  ,  within  and  ad¬ 
jacent  to  the  contact  zone.  Within  the  contact  zone  (x/b  <_  1.00)  the 
classical  pressure  distribution  expected  beneath  a  conical  indenter  is 
obtained  from  the  elastic  finite  element  analysis  for  the  prescribed 
displacement  loading.  As  can  be  seen  from  the  figure,  this  corresponds 
very  well  to  that  obtained  from  the  closed-form  Sneddon  solution  for 
elastic  conical  indentation.  This  verifies  the  displacement  loading 
condition  assumed  in  this  investigation,  where  only  vertical  dis¬ 
placements  of  surface  nodal  points  are  specified.  In  addition,  the 
load  resultant  predicted  from  the  Sneddon  solution  (Appendix  A), 

F  =  3Q04N,  compares  well  with  that  determined  independently  by 
s 

summation  of  the  vertical  nodal  forces  on  the  contact  surface  of 
0  <  x/b  <_  1.00,  F®  -  2943N.  This  factor  further  substantiates  the 
displacement  method  of  load  application. 


*  Sneddon  Analysis 

■  Elastic  Finite  Element  Analysis 

■  Elastic/Plastic  Finite  Element  Analysis 


Figure  21.  Non-dinensionalized  Pressure  Distribution 

on  and  adjacent  to  the  Contact 


4.3.2  Elastic/Plastic  Conical  Indentation 

In  generating  elastic/plastic  stress  states  beneath  the  conical 
indenter,  the  model  was  allowed  to  yield  according  to  the  Huber- 
Mises  yield  criterion  at  nodal  points  beneath  the  indenter  with  sub¬ 
sequent  flow  being  governed  by  the  Prandtl-Reuss  flow  rule.  The 
elastic/perfectly  plastic  stress-strain  curve  was  assumed  for  soda- 
lime  glass.  During  each  iteration  loop,  both  the  elastic  and  elastic/ 
plastic  portions  of  the  stiffness  matrix  were  updated.  As  stated 
earlier,  the  vertical  nodal  point  displacement  along  the  contact  sur¬ 
face  was  used  as  load  input  for  an  assumed  contact  radius  of  0.30  mm. 

The  elastic  Sneddon  solution  was  used  to  predict  vertical  displacements, 
thus  raising  an  obvious  question  regarding  the  use  of  an  elastic 
solution  for  the  prediction  of  elastic/plastic  displacements  beneath 
the  conical  indenter. 

The  applicability  of  the  elastic  Sneddon  solution  for  the  pre¬ 
diction  of  vertical  displacements  on  the  contact  surface  was  judged 
by  observation  of  mesh  deformation  for  the  elastic/plastic  case,  and 
the  extent  and  character  of  the  plastic  region.  The  mesh  deformation 
for  the  elastic/plastic  conical  indentation  is  shown  in  Figure  22. 

The  corresponding  extent  and  character  of  the  plastically  deformed 
region  and  the  stress  trajectories  beneath  the  indenter  can  be  seen  in 
Figure  23.  These  figures  indicate  a  continuity  of  surface  deformation 
both  within  and  adjacent  to  the  contact  surface,  and  also  plastic 
region  characteristics  similar  to  those  measured  experimentally  and 
predicted  analytically.  The  shape  of  the  plastic  zone  obtained  by 


Figure  22.  Elastic/Plastic  Mesh  Deformation. 
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Shaded  Area:  Plastic  Zone 
< - *  :  (Tensile  Stress) 

>*—  :  o ^  (Compressive  Stress) 

<-c-»  :  a2  (Tensile  Stress)  1 

^  1  (a^  is  Out-of- 

„2  (Compressive  Stress) •  Plane  Hoop  Stress) 

Figure  23.  Plastic  Zone  and  Elastic/Plastic  Stress 
Traj  ectories . 


the  elastic/plastic  finite  element  analysis  is  nearly  hemispherical 
and  is  centered  directly  beneath  the  indenter.  Its  radius  c  is 
approximately  0.35  mm,  resulting  in  a  ratio  of  c/b  equal  to  1.16. 

This  c/b  ratio  compares  favorably  with  a  value  of  1.159  explicitly 
obtained  by  Johnson  [11]  (APPENDIX  B)  as  a  function  of  material 
properties,  with  the  center  of  the  plastic  region  at  the  undeformed 
origin.  These  factors  were  therefore  considered  as  substantiation 
of  the  assumed  displacement  loading  as  obtained  from  the  Sneddon 
solution.  Observation  of  Figure  23  illustrates  the  difference  in  the 
principal  stress  trajectories  of  the  elastic  and  plastic  regions.  In 
the  elastic  region,  the  in-plane  hoop  stress  is  the  maximum  principal 
stress,  a^,  whereas  in  the  plastic  region,  the  out-of-plane  hoop  stress 
is  mainly  the  maximum  principal  stress,  cr^. 

In  figures  17  to  20,  the  elastic/plastic  stress  fields  along 
the  surface  of  the  half-space  are  compared  with  the  elastic  solutions 
obtained  from  the  finite  element  analysis  and  the  Sneddon  solution. 

The  horizontal  stress,  a  ,  along  the  surface  shows  the  characteristic 

XX 

tensile  peak  at  the  elastic/plastic  boundary  as  indicated  in  Figure  17. 
Within  the  contact  zone,  this  stress  is  again  compressive,  reaching  a 
compressive  peak  directly  beneath  the  indenter.  Outside  the  region 
of  contact,  this  stress  rapidly  approaches  the  value  obtained  in  the 
elastic  analysis  at  the  elastic/plastic  boundary.  This  tensile 
peak  is  usually  associated  with  ring  cracking,  as  suggested  in  the 
previous  section.  The  out-of-plane  hoop  stress,  ®zz*  along  the  sur¬ 
face  again  is  everywhere  compressive,  with  the  peak  compression  being 
found  directly  beneath  the  indenter,  as  indicated  in  Figure  18. 
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This  stress  does  not  rapidly  approach  elastic  results  with  increasing 
distance  from  the  contact  region  as  was  seen  in  the  case  of  the  hori¬ 
zontal  stress,  a 

xx 

Directly  beneath  the  indenter,  the  horizontal  stress,  a  ,  changes 
from  high  compressive  values  near  the  surface  to  small  tensile  values 
at  the  elastic/plastic  boundary,  as  indicated  in  Figure  19.  This 
stress  rapidly  approaches  levels  obtained  in  the  elastic  analysis 
near  the  elastic/plastic  boundary.  Tensile  values  of  this  stress 
result  in  median  cracking  as  suggested  in  the  previous  section.  The 
subsurface  vertical  stress  component,  a  ,  i-s  again  everywhere  com¬ 
pressive;  reaching  a  peak  magnitude  directly  beneath  the  indenter. 

It  decreases  to  zero  with  increasing  distance  from  the  surface  more 
rapidly  than  the  corresponding  stress  in  the  elastic  analysis.  This 
is  indicated  in  Figure  20. 

In  Figure  21,  the  elastic/plastic  pressure  distribution,  o  , 
within  and  adjacent  to  the  contact  zone  ic  compared  with  the  elastic 
solutions  obtained  from  the  finite  element  analysis  and  Sneddon 
solution.  It  is  obvious  that  high  pressure  exists  within  the 
region  of  contact  and  that  this  pressure  is,  of  course,  suddenly 
reduced  beyond  the  contact  surface.  The  elastic/plastic  load 
resultant  generated  from  summation  of  the  vertical  nodal  forces  on 
the  contact  surface  of  0  <  x  <  b  is  obtained  as  F^  =  2498N.  This 
is  15  to  17  percent  lower  than  the  elastic  load  resultants  obtained 
from  the  Sneddon  solution  and  the  elastic  finite  element  analysis. 

This  lower  load  resultant  value  was  expected,  since  less  pressure 
was  needed  to  deform  the  yielded  region. 


CHAPTER  V 


RESIDUAL  SOLUTIONS  FOR  ELASTIC/PLASTIC  CONICAL  INDENTATION 

5. 1  Introduction 

In  this  chapter,  residual  stress,  displacement  and  strain  fields 
resulting  from  the  elastic/plastic  conical  indentation  are  generated. 
Development  of  the  residual  solution  required  that  a  load-unload  cycle 
consisting  of  two  separate  loading  stages  (stage  II  and  III)  be  formu¬ 
lated.  During  the  loading  stage  of  the  cycle  (stage  II),  equivalent 
nodal  forces  were  applied  as  load  input  to  surface  nodes  within  the 
contact  area.  These  equivalent  nodal  loads  were  obtained  from  the  nodal 
force  output  obtained  in  stage  I  displacement  loading  and  produce  the 
surface  deformation  and  resultant  stress,  displacement  and  strain 
distribution  beneath  the  indenter.  A  procedure  developed  to  convert  the 
nodal  force  output  obtained  in  stage  I  to  the  equivalent  nodal  load  is 
discussed  in  greater  detail  in  the  following  sections  of  this  chapter. 

In  the  unloading  stage  of  the  cycle  (stage  III),  the  equivalent  nodal 
loads  were  reduced  to  zero,  resulting  in  the  required  residual  solutions. 


5.2  Method  of  Analysis 


The  procedure  used  to  obtain  the  equivalent  nodal  loads  in  the 


contact  area  from  the  nodal  force  output  obtained  in  stage  I  is  dis¬ 
cussed  in  this  section.  These  equivalent  nodil  loads  are,  as  mentioned 
previously,  used  during  the  loading  stage  of  the  cycle  (stage  II). 


In  order  to  improve  stress  and  displacement  fields  beneath  the 
indenter  for  stage  II  loading,  an  iterative  scheme  was  developed  through 
which  the  output  nodal  forces  from  stage  I  loading  were  continuallv  im¬ 
proved  by  incrementing  each  nodal  load  magnitude  until  the  error  in  re¬ 
sultant  displacement  and  stress  fields  and  nodal  force  resultants  be¬ 
neath  the  indenter  was  less  than  3  percent.  With  the  error  reduced  to 
less  than  3  percent,  the  nodal  force  loading  was  considered  to  be 
equivalent  to  the  nodal  displacement  loading.  Figure  24  compares  the 
nodal  forces  resulting  from  the  displacement  loading  to  the  equivalent 
nodal  loads  obtained  through  the  iterative  scheme. 

Following  t  ie  equivalent  nodal  force  loading  in  the  loading  phase 
of  the  cycle,  the  equivalent  nodal  loads  were  reduced  to  zero  yielding 
residual  displacements,  strains  and  stresses  on  the  contact  surface 
and  beneath  the  indentation. 

5.3  Results  and  Discussion 

A  geometric  comparison  of  the  indentation  profile  obtained  in 

the  loading  phase  of  the  cycle  and  the  residual  surface  crater  remaining 

after  unloading  is  shown  in  Figure  25.  From  Figure  25,  the  ratio  of 

residual  crater  depth  to  maximum  indentation  depth,  v  /v  ,  is 

easily  obtained  as  0.41.  This  ratio  can  be  compared  to  that  analytically 

predicted  from  the  equation  developed  by  Lawn  and  Howes  [46]  (APPEND LX 
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B,  equation  (B-8)).  The  analytical  ratio  obtained  is  v  /v  =  0.70 

oiflx  max 

for  a  hardness  of  soda-lime  glass  of  5.6  GPa  and  the  assumption 

y  =  y  ,  where  and  y„  are  elastic  and  plastic  hysteresis  geometrical 
EH  E  H 

parameters,  respectively  (APPENDIX  B,  equation  (B-9)).  The  disagreement 


Nodal  Forces  resulted  by  Indentation  (Output  of 
Stages  I) 

Equivalent  Nodal  Loads  (Input  for  Stage  II) 


24.  Comparison  of  Nodal  Forces  Resulting  from 
Indentation  and  Equivalent  Nodal  Loads. 


ini; 


o  =  Indentation 
o  =  Residual  Crator 


Figure  25.  Geometrical  Comparison  of  Indentation 
Profile  and  Residual  Crater. 
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between  the  finite  element  and  analytical  results  may  be  caused  by  the 
assumption  y_  =  yu.  By  using  the  value  of  v  /v  =  0.41  predicted 
by  the  finite  element  solution  and  equation  (B-8) ,  a  geometrical  para¬ 
meter  ratio  Yjj/Yjj  of  1.26  is  obtained  resulting  in  a  geometrical  para¬ 
meter  for  plastic  hysteresis  y^  of  1.25.  This  result  is  based  on  an 
assumed  elastic  geometrical  parameter  y  =  ir/2  defined  by  Sneddon  [5,6]. 

It  should  be  noted  that  there  is  a  pile-up  on  the  free  surface 
adjacent  to  the  contact  radius,  as  shown  in  Figure  25,  and  this  pile-up 
is  the  result  of  elastic  recovery  of  the  surface  elements  adjacent  to 
the  contact  area.  Figure  26. 

Non-dimensionalized  horizontal  and  out-of-plane  hoop  residual 

stresses  along  the  indentation  surface  are  shown  in  Figure  27. 

Stresses  are  shown  as  a  function  of  non-dimensionalized  distance  from 

the  center  of  indentation.  Both  the  horizontal  and  out-of-plane  hoop 
R  R 

stresses,  a  and  o  ,  are  tensile  for  the  regions  0.5  <  x/b  <1.5 

xa  Z  Z  " 

and  x/b>  0.6,  respectively,  with  tensile  peaks  occurring  at  the 
elastic/plastic  boundary  x/b  =  1.3.  These  tensile  residual  stresses 
may  be  responsible  for  the  propagation  of  radial  cracking  during  the 
unloading  phase  of  the  cycle,  a  topic  of  discussion  in  a  subsequent 
chapter. 

Non-dimensionalized  residual  stresses  directly  beneath  the  indenter 


are  shown  in  Figure  28.  The  horizontal  and  out-of-plane  hoop  stress, 

R  R  R 

a  and  a  ,  as  well  as  the  vertical  stress  component,  <j  ,  are 


tensile  for  y/b  >  1.5.  Tensile  values  of  the  out-of-plane  hoop  residual 


stress  could  result  in  median  crack  extension  during  the  unloading 


phase.  The  vertical  residual  tensile  stresses  could  promote  lateral 


•a*l 


[im 


A  -  a^x/Etang 
o  *  o^/Etang 


Figure  27. 


Non-dimensionalized  Residual  Surface 
Stresses  along  the  x-axis. 
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crack  propagation  upon  unloading.  Both  of  these  possibilities 
will  be  explored  in  the  following  chapter. 

The  residual  plastic  zone  and  stress  trajectories  are  shown  in 
Figure  29.  Comparison  with  Figure  23  reveals  that  the  size  of  the 
plastic  zone  has  increased  approximately  20  percent  during  the  un¬ 
loading  phase  due  to  elastic  recovery  in  the  plastic  region.  Of 
major  interest  are  the  stress  trajectories  in  the  elastic  region 
where  fracture  can  be  extended  during  the  unloading  phase.  In 
this  region  the  stress  trajectories  indicate  the  possible  propagation 
of  the  lateral  cracking  driven  by  the  in-plane  hoop  stress,  o^, 
from  points  near  the  elastic/plastic  boundary.  Near  the  surface, 

the  trajectories  indicate  that  the  out-of-plane  hoop  residual 
R  R 

stress,  a  *  a.,  would  result  in  the  radial  crack  extension.  For 

ZZ  X 

subsurface  regions,  the  trajectories  indicate  the  possibility  of 

the  median  crack  propagation  by  the  out-of-plane  hoop  residual  stress, 
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Figure  29.  Plastic  Zone  and  Residual  Stress  Trajectories. 


CHAPTER  VI 


ANALYSIS  OF  FRACTURE  MECHANISMS  INDUCED  BY  INDENTATION 

6.1  Introduction 

Analytical  and  experimental  investigations  considering  the 
mechanics  of  crack  initiation  and  propagation  beneath  spherical, 
conical,  Vickers  pyramidal,  and  point  indentations  in  a  variety  of 
materials  have  been  the  topic  of  many  recent  investigations  [29-48]. 

In  this  chapter,  fracture  mechanisms  beneath  conical  indentation  in 
soda-lime  glass  are  further  described  by  applying  the  numerical  results 
obtained  in  the  preceding  chapters  to  approximate  fracture  mechanics 
theories.  Conical  indentation  was  chosen  for  analysis  as  it  closely 
simulates  loading  conditions  found  in  a  large  variety  of  machining 
processes. 

6.2  Failure  Modes  beneath  the  Indenter 
There  are  four  major  crack  systers  introduced  as  indentation 

fracture  mechanisms.  These  crack  systems  can  be  classified  by  the 
cracks  initiated  and  propagated  near  the  surface  and  at  subsurface 
locations.  It  has  also  been  experimentally  and  analytically  show  that 
these  crack  systems  are  initiated  on  or  near  the  elastic/plastic  boundary 
with  subsequent  extension  into  the  elastic  region.  The  schematic 
diagram  of  these  crack  systems  is  shown  in  Figure  30. 
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6.2.1  Crack  Systems  near  the  Surface 

Near  the  surface,  two  major  crack  systems  are  usually  found  in 
axisymmetric  indentation.  They  are  surface  ring  crack  system  which  is 
initiated  in  the  loading  phase  of  the  cycle  [35,36],  and  radial  crack 
system  which  is  usually  initiated  and  propagated  during  the  unloading 
phase  of  the  cycle  [33,35,36].  It  is  known  that  surface  ring  crack 
system  is  often  propagated  as  a  Hertzian  cone  crack  [31,35,37,40]  in 
the  case  of  spherical  indentation.  In  the  case  of  Vickers  pyramidal 
indentation,  radial  crack  system  is  often  initiated  in  the  loading  phase 
and  propagated  during  the  unloading  phase  of  the  cycle,  since  the  edges 
of  pyramidal  indent er  create  highly  singular  regions  in  the  specimen 
[33,36,41,42]. 

6.2.2  Subsurface  Crack  System 

Two  major  subsurface  crack  systems  are  also  found.  They  are  me¬ 
dian  crack  system,  which  is  initiated  and  propagated  during  the  load¬ 
ing  phase  and  further  propagated  during  the  unloading  phase  of  the 
cycle  [33,35,36,38],  and  lateral  crack  system,  which  is  usually  ini¬ 
tiated  and  propagated  during  the  unloading  phase  of  the  cycle.  It 
is  known  that  the  shape  of  the  median  cracking  is  penny- like  or  half¬ 
penny-like  [33,34,36,38,39,43]  beneath  the  indentation.  Lateral  crack 
system  is  often  initiated  with  the  surface  ring  crack  and  propagated 
as  a  Hertzian  cone  crack  in  the  case  of  spherical  indentation  [31,35, 
37,40]  as  described  in  the  last  section.  In  the  case  of  Vickers  pyra¬ 
midal  indentation,  lateral  crack  system  is  initiated  on  the  elastic/ 
plastic  boundary  and  propagated  to  the  surface  [33,35,37,38.41]. 
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6. 3  Fracture  Initiation  and  Propagation 

In  the  discussion  of  fracture  initiation  and  propagation  during 
elastic/plastic  conical  indentation,  only  the  opening  mode  (K^-mode) 
was  considered.  To  correlate  the  stresses  numerically  obtained  by 
the  finite  element  method  with  opening  mode  fracture,  the  tensile 
peak  stresses  at  the  surface  and  at  subsurface  locations  must  be 
considered,  especially  near  the  elastic/plastic  boundary. 

6.3.1  Identification  of  the  Tensile  Peak  Stresses 

The  necessary  tensile  stress  fields  and  tensile  stress  trajectories 
for  loading  and  unloading  of  a  conical  indenter  were  numerically 
determined  in  the  preceding  chapters.  By  observation  of  Figures  17 
and  27,  the  maximum  surface  tensile  stresses  are  found  at  the  elastic/ 
plastic  boundary.  In  the  loading  phase,  the  maximum  value  of  the 
horizontal  surface  stress,  >  corresponding  to  surface  ring 
cracking,  is  obtained  as  2.2  GPa.  During  the  unloading  phase,  the 

D 

maximum  values  of  the  horizontal  surface  residual  stress,  a  , 

xx 

corresponding  to  the  surface  ring  cracking,  and  the  surface  out-of- 

£ 

plane  residual  hoop  stress,  a  ,  corresponding  to  radial  cracking, 

z  z 

are  obtained  as  1.6  GPa  and  1.8  GPa,  respectively.  By  observation  of 

Figures  19  and  28,  the  maximum  subsurface  tensile  stresses  along 

the  y-axis  are  also  found  at  the  elastic/plastic  boundary.  In  the 

loading  phase,  the  maximum  value  of  the  subsurface  horizontal  stress, 

a  ,  which  is  equal  to  the  subsurface  out-of-plane  hoop  stress,  n  , 

XX  z  z 

corresponding  co  the  median  cracking,  is  obcained  as  0.48  GPa.  During 
the  unloading  phase  of  the  cycle,  the  maximum  values  of  the  subsurface 


66 


horizontal  residual  stress,  a  =  a  ,  corresponding  to  the  median 

XX  2  Z 

cracking,  and  the  subsurface  vertical  residual  stress,  ctR  ,  corres- 

yy 

ponding  to  the  lhterial  cracking,  are  obtained  as  0.42  GPa  and  0.24  GPa, 
respectively.  Those  tensile  peak  stresses  are  shown  in  Table  1  with 
the  corresponding  crack  systems  which  they  initiate. 

6.3.2  Critical  Flaw  Size 

The  peak  tensile  stresses  obtained  in  the  previous  section  can 
be  combined  with  an  appropriate  fracture  mechanics  theory  to  predict 
fracture  initiation  and  propagation  beneath  the  indenter.  Crack 
initiation  is  here  predicted  by  the  Griffith  flaw  hypothesis  [29,30]: 

a  >.(2Er/ircf)^  ,  (6,1) 

where  a  is  the  critical  stress  causing  stable  fracture,  E  is  Young's 
modulus,  T  is  the  crack  surface  energy  and  cf  is  the  flaw  size.  By 
solving  equation  (6.1)  for  c^,  the  following  expression  is  obtained 
for  the  critical  flaw  size  corresponding  to  the  tensile  peak  stress: 

cf  >  (2Er/7ro2).  (6.2) 

On  the  surface,  the  critical  flaw  sizes  corresponding  to  the 
horizontal  tensile  peak  stress,  2.2  GPa,  and  the  horizontal  residual 
peak  stress,  1.6  GPa,  are  obtained  from  equation  (6.2)  as  0.025  urn 
and  0.048  pm,  respectively.  These  critical  flaw  sizes  will  result  in 
initiation  and  propagation  of  the  surface  ring  cracking.  Also  on  the 
surface,  the  critical  flaw  size  corresponding  to  the  out-of-plane  hoop 
residual  tensile  peak  stress,  1.8  GPa,  can  be  calculated  as  0.039  pm 
which  will  result  in  initiation  and  propagation  of  the  median  cracking. 


Table  1.  Tensile  Peak  Stresses  and  the  Corresponding  Crack  Systems 
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At  subsurface  locations  along  the  y-axis,  the  critical  flaw  sizes 
corresponding  to  the  horizontal  or  out-of-plane  hoop  tensile  peak  stress, 
0.48  GPa,  and  the  horizontal  or  out-of-plane  hoop  residual  tensile  peak 
stress,  0.42  GPa,  are  obtained  from  equation  (6.2)  as  0.54  um  and  0.7  uni, 
respectively.  These  critical  flaw  sizes  will  result  in  initiation  and 
propagation  of  median  crack  system.  Also  at  subsurface  locations  along 
the  y-axis,  the  critical  flaw  size  corresponding  to  the  vertical  residual 
tensile  peak  stress,  0.24  GPa,  can  be  obtained  as  2.17  um,  which  will 
result  in  initiation  and  propagation  of  lateral  crack  system. 

The  critical  flaw  sizes  have  been  discussed  in  recent  papers. 

Lawn  and  Evans  [36]  have  given  the  critical  condition  for  the  initiation 
of  subsurface  flaws  as: 

cf  =  44.2  (KIC/H)2  .  (6.3) 

Hagan  [43,44]  has  developed  a  similar  expression  for  the  largest  flaw 
size  which  could  be  nucleated  by  the  intersecting  flow  lines  beneath 
the  indenter: 

cf  =  2.95  (KIC/H)2  .  (6.4) 

In  these  equations,  and  H  are  the  critical  stress  intensity  factor 

-3/2 

and  hardness,  respectively.  For  values  of  =  0.7  IINm  and 

H  =  5.6  GPa  for  soda-lime  glass,  equations  (6.3)  and  (6.4)  predicted 
critical  flaw  sizes  of  0.7  um  and  0.46  um,  respectively.  The  void  size 
for  soda-lime  glass  has  been  experimentally  measured  by  Hagan  [44] 
as  0.6  um.  These  flaw  sizes  were  correlated  to  the  median  cracking. 
Swain  and  Hagan  [35]  and  Chiang,  Marshall  and  Evans  [47]  correlated 


the  crack  systems  and  the  threshold  flaw  sizes.  The  threshold  flaw 
sizes  for  median  and  radial  crack  systems  were  estimated  as  1.0  ym 
and  that  for  lateral  crack  system  was  estimated  as  3.0  am.  All  of  those 
critical  flaw  sizes  are  compared  in  Table  2.  The  critical  flaw  size 
obtained  by  the  finite  element  analysis  with  the  Griffith  energy 
hypothesis  agrees  very  well  with  experimental  and  analytical  values 
for  median  crack  system.  However  the  critical  flaw  sizes  obtained  by 
Swain  et  al.  and  Chiang  et  al.  are  clearly  overestimated  as  compared 
with  other  analytical,  experimental  and  numerical  values  for  crack 
initiation, 

6.3.3  Fracture  Initiation  and  Propagation  Zones 

As  described  previously,  the  critical  flaw  sizes  required  for 
the  initiation  and  propagation  of  the  various  crack  systems  are 
obtained  from  the  Griffith  energy  hypothesis  with  the  stresses  generated 
by  the  finite  element  method.  It  was,  then,  found  that  the  critical 
flaw  sizes  predicted  by  Swain  et  al.  and  Chian  et  al.  were  overesti¬ 
mated  when  compared  with  the  critical  flaw  sizes  obtained  in  this 
investigation. 

However,  by  use  of  these  overestimated  flaw  sizes,  the  critical 
stresses  can  be  generated  by  equation  (6.1)  to  predict  probable 
crack  initiation  and  propagation  zones.  The  values  of  the  critical 
stress  corresponding  to  initiation  and  propagation  of  the  median  and 
radial  crack  systems  is  obtained  as  0.35  GPa  from  equation  (6.1)  with 
the  value  of  c^  =  1.0  um.  Also  the  value  of  the  critical  stress 
corresponding  to  lateral  crack  system  can  be  obtained  as  0.2  GPa  from 
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equation  (6.1)  with  the  value  of  =  3.0  um.  These  two  critical 
stress  values  can  be  compared  with  the  stresses  generated  by  the 
finite  element  analysis.  In  Figure  31,  the  assumed  crack  initiation 
and  propagation  zones  corresponding  to  a  >_  0.35  GPa  in  the  loading 
phase  are  plotted.  Those  can  be  thought  of  as  the  assumed  surface 
ring  crack  zone  and  the  assumed  median  crack  zone.  In  Figure  32, 
the  assumed  crack  initiation  and  propagation  zones  corresponding 
to  o_>0.35  GPa  and  >_  0.2  GPa  during  the  unloading  phase  of  the  cycle 
are  plotted.  Those  can  be  thought  of  as  the  potential  median  crack, 
radial  crack  and  lateral  crack  zones. 
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axx  >0.35  GPa  (Potential  Surface-Ring  Crack  Zone) 


0.35  GPa  (Potential  Median  Crack  Zone) 


Figure  31.  Crack  Initiation  and  Propagation  Zones 
in  the  Loading  Phase 
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a  \  3_  0.35  GPa  (Potential  Radial  Crack 

zz  Zone) 


(Potential  Lateral  Crack 
Zone) 


Figure  32.  Crack  Initiation  and  Propagation  Zones  during 
the  Unloading  Phase. 


CHAPTER  VII 


CONCLUSIONS  AND  RECOMMENDATIONS 


7. 1  Conclusions 

The  following  conclusions  were  reached  in  this  investigation: 

1.  The  elastic/plastic  finite  element  code  used  in  this  investi¬ 
gation  was  verified  by  comparing  nodal  displacements,  forces 
and  plastic  zone  extension  obtained  for  selected  plane- 
stress  problems,  the  plates  with  a  circular  penetration,  with 
results  numerically  obtained  by  Yamada  [Personal  Note]  and 
Marcal  and  King  [18]  in  independent  investigations. 

2.  The  solution  technique  was  further  verified  by  comparison 
of  elastic  displacements,  strains  and  stresses  generated 
for  conical  indentation  of  an  elastic  half-space  to  the 
closed-form  results  obtained  from  the  Sneddon  solution. 

3.  The  elastic/perfectly  plastic  finite  element  model  developed 
for  conical  indentation  of  a  half-space  was  successfully 
utilized  to  determine  the  size  and  characteristics  of  the 
plastic  zone  beneath  the  indenter  during  a  loading-unloading 
cycle.  The  size  of  the  plastic  zone  was  shown  to  agree 
well  with  results  obtained  analytically  by  Johnson  [11]. 

The  shape  of  the  plastic  zone  was  shown  to  be  hemispherical 
which  conformed  to  results  obtained  in  previous  experimental 


and  analytical  investigations.  The  origin  of  the  hemis¬ 
pherical  plastic  zone  was  found  to  be  on  the  deformed 
surface  directly  beneath  the  indenter,  not  on  the  undeformed 
surface. 

The  displacement,  strain  and  stress  fields  were  successfully 
generated  by  the  elastic/plastic  finite  element  model  for  a 
loading-unloading  cycle.  Resulting  loading  and  residual 
stress  fields  were  examined  to  determine  the  location  and 
magnitude  of  tensile  peak  stresses  which  would  result  in 
the  initiation  and  subsequent  propagation  of  the  ring,  radial, 
median  and  lateral  crack  systems.  Tensile  peak  stresses 
were  characteristically  found  on  the  elastic/plastic  boundary 
near  the  surface  and  at  subsurface  locations  directly 
beneath  the  indentation.  The  out-of-plane  hoop  stress,  which 
controlled  both  radial  and  median  cracking,  was  found  to 
change  sign  near  the  surface  from  compression  during  the 
loading  phase  to  tension  during  the  unloading  phase  of  the 
cycle.  This  indicates  that  radial  crack  propagation 
and  the  additional  median  crack  propagation  takes  place 
during  the  unloading  phase,  a  fact  experimentally  verified. 

By  use  of  the  Griffith  energy  hypothesis  with  the  generated 
tensile  peak  stresses,  the  critical  flaw  sizes,  resulting  in 
initiation  of  the  crack  systems,  were  obtained.  The  critical 
flaw  sizes  obtained  in  this  investigation  agreed  well  with 


critical  flaw,  sizes  predicted  analytically  and  experimentally 
by  Lawn  et  al.  [36],  Hagan  [43,44],  Strain  et  al.  [35]  and 
Chiang  et  al.  [47]  for  each  crack  system. 

6.  Fracture  initiation  and  propagation  zones  for  the  loading 

phase  and  unloading  phase  of  the  cycle  were  generated  through 
comparison  of  the  tensile  peak  stresses  obtained  in  this 
investigation  to  critical  stresses  obtained  from  the  Griffith 
energy  hypothesis  and  maximum  allowable  flaw  sizes.  These 
fracture  initiation  and  propagation  zones  agreed  well  with 
actual  initiation  and  propagation  zones  experimentally 
observed. 

7.2  Recommendat ions 

The  following  subjects  are  recommended  as  topics  for  further 
research: 

1.  The  elastic/plastic  finite  element  model  should  be  expanded 
to  consider  elastic/plastic  strain  hardening  in  order  to 
analyze  more  realistic  indentation  models. 

2.  Investigation  should  be  expanded  to  consider  more  brittle, 
lower  yield  strength/Young's  modulus  ratio  models  such  as 
MgO  and  ZnS  in  order  to  determine  the  effect  of  these  para¬ 
meters  on  the  extent  and  character  of  the  plastic  zone  and 
resultant  crack  systems. 

3.  Shear  peaks  as  well  as  tensile  peaks  should  be  investigated 
as  they  affect  fracture  initiation  through  dislocation 
activity  and  crack  propagation  by  Mode  II  loading. 
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4.  The  elastic/plastic  finite  element  model  should  be  expanded 
to  consider  non-axisymmetric  indentation  such  as  Vickers 
pyramidal  indentation  and  obtain  characteristics  of  the 
singularities  created  by  the  edges  of  the  indenter. 

5.  The  elastic/plastic  finite  element  technique  should  be 
combined  with  statistical  analysis,  such  as  "The  Monte 
Carlo  method','  in  order  to  closely  model  elastic/plastic 


indentation  fracture. 


APPENDIX  A 


THE  ELASTIC  STREAS  FIELDS  FOR  POINT,  SPHERICAL  AND  CONICAL  INDENTATION 

1.  Point  Indentation  (Boussinesque  and  Michell  Solution) 

The  well  known  non-dimens ionalized  stress  field  equations  for 
the  point  indentation  problem  [2]  are  given  for  the  coordinate 
system  shown  in  Figure  33  as: 

V  =  2(i T  {‘  (3X  +  4u)cos6  +  y} 

°  a  V  2„  2  8 

7?  ■  COS  9  sec  2  • 

^  IT  (cos9  "  sin20)  sec2  f  * 

F9  =  40^0  Sine  COS0  sec2  1  ’ 

Tr<f>  =  Te*  "  °’ 
where 

X  =  _ I _ 

(1+v) (l-2v)  ’ 

E 

U  =  2 (1+v) 


F  *  applied  load 


2.  Spherical  Indentation  (Hertz  and  Huber  Solution) 

The  non-dimensionalized  Hertzian  stress  field  which  was  extended 
to  a  general  three-diinensional  form  by  Huber  [4]  is  expressed  by  four 
sets  of  equations  depending  on  the  coordinates  of  the  point  in  the 
half-space.  Figure  34. 


.  For  points  along  the  y-axis  directly  beneath  the  indenter 
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For  points  on  the  surface  of  the  half-space  beneath  the  contact  sur¬ 


face  of  the  indenter  (y=0,  R  <_  a) : 


\  l-2v  a2  M  ,1/2  „2,3,  1.5  ,  2  „2 

r  =  —  3  [1  -  (-  /a  -R  )  ]  -  —  /a  -R  , 

°  y=0  R 


=  -1.5 


/a2-R2 


¥4  [1  -  (t  /a2-R2)3]  -  /a2-R2 


AD-A127  529  AN  ANALVSIS  OF  CONICAL  INDENTATION  OF  AN 

ELASTIC/PERFECTLV  PLASTIC  HALF-S.  .  <U>  UENNSVL VANIA 
STATE  UNIV  UN  I VERS I TV  PARK  APPLIED  RESEARCH  LAB.  . 
UNCLASSIFIED  T  KATO  26  OCT  82  ARL/PSU/TN-82-230  F/G  28/11 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS-1963-A 


82 


and 


8Ea 

3D(l-v2) 


3.  Conical  Indentation  (Sneddon  Solution) 

The  coordinate  system  for  Sneddon's  analysis  is  given  in  Figure 
35.  Before  defining  expressions  for  the  stress  distributions  in 
the  half-space,  the  displacement  functions,  mean  pressure  and  pressure 
distribution  at  the  surface  of  the  half-space  will  be  defined.  The 
vertical  displacement  at  the  surface  of  the  half-space  is  assumed  as: 
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The  non-dimensionalized  stresses  within  the  half  space  are  given  by 
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K2  =  p  [ Q  cos^-^>  “  '  -  1  , 


/P2  +  C2 


K0  =  2  tpKl  ‘  ai  "  p  (1'Q  Sin*)]  ’ 


and 

V  X 

C  *  £  »  p=b’  e  =  b  tan  8  , 

2  2  .4  .  2  .  2  n2  2 

t  =  H5  ,  Q  -  (p  +  5  -1)  +4^  , 

<|>  =  arctan  (~)  »  $  ■  y  arc  tan 
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(p2+52-l)| 


4.  Comparison  of  the  three  elastic  indentation  stress  fields. 

To  compare  the  three  elastic  stress  fields  which  are  given  in  the 
preceding  sections  of  APPENDIX  A  for  the  same  condition.  Young's 
modulus  E,  Poisson’s  ratio  v,  and  the  resultant  force  F  must  be  fixed. 
For  the  finite  element  analysis  of  conical  indentation,  the  following 
four  values  were  defined: 

E  -  49000  N/ram2  v  -  0.22 

Contact  radius  b  =  0.3  mm  cone  angle  a  =  136°  (8  =  22°) 

From  these  four  values  the  resultant  force  F  can  be  obtained  from 


equation  (A-8): 


2(l-v) 


If  we  equate  the  contact  radius  a  for  Huber's  with  b  for  Sneddon's 
then  the  spherical  diameter  of  Huber  solution  can  be  defined  as: 

D  =  — ■  =  14.00  (mm)  . 

3F(l-v  ) 

Huber's  stress  fields  and  Sneddon's  stress  fields  were  non-dimens ionalied 
2  2 

by  the  factor  tra  /F  =  rb  /F;  however,  Michell's  stress  field  was  non- 

2 

dimensionalized  by  the  factor  irr  /F  (r  is  variable).  To  compare  Michell's 
stress  field  with  other  two  stress  fields  results  for  the  Michell  solu- 
tion  (A-l)  are  multiplied  by  p’  *  F/rrr  and  divided  by  Pp=pm=F/Tra  =F/ub  . 
The  comparison  of  the  three  elastic  non-dimens ionalized  stress  fields  in 
terms  of  the  angle  0  are  plotted  in  Figures  36  through  39  for  the 
coordinates: 

2  2  2  2 
x  +  y  *  a  =  r  and  z  =  0, 

where  x  -  rsin0  ,  y  =  rcos0  . 

To  see  the  effect  of  the  cone  angle  a  in  Sneddon  solution,  the 

cone  angle  a  and  the  contact  radius  b  are  assumed  as  variables  while 

E,  v,  F,  a  and  D  are  remained  as  constants.  The  relationship  between 

a, 8,  b  are  shown  in  Table  3.  With  the  constant  non-dimensionalized 

factor  l/pQ>  the  change  of  Sneddon's  stress  field  in  terms  of  the  angle  9 

are  plotted  in  Figures  40  through  43  under  the  condition: 


ssss£HbbSh 


Angle  e( 8 ) 

Michell  Solution  (Point  Indentation) 
Huber  Solution  (Spherical  Indentation) 
Sneddon  Solution  (Conical  Indentation) 


Figure  35.  Comparison  of  Elastic  Point,  Spherical  and 
Conical  Indentation  Stress  Fields  (I): 
Non-dimensionalized  Stress,  a  ... 


Angle  0(8) 

M  :  Michell  Solution  (Point  Indentation) 

H  :  Huber  Solution  (Spherical  Indentation) 
S  :  Sneddon  Solution  (Conical  Indentation) 


Figure  38.  Comparison  of  Elastic  Point,  Spherical 

and  Conical  Indentation  Stress  Fields  (III) 
Non- dimens ionalized  Stress,  a _ 
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Figure  39. 


Comparison  of  Elastic  Point,  Spherical 

and  Conical  Indentation  Stress  Fields  (IV) 

Non- dimens ionalized  Stress,  a  . 
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Angle  9 ( 8 ) 

A:  Sneddon  Solution  (Conical  Indentation) :  a  =  160° 

B:  Sneddon  Solution  (Conical  Indentation):  a  =  120° 

C:  Sneddon  Solution  (Conical  Indentation):  a  =  80° 

D:  Sneddon  Solution  (Conical  Indentation) :  a  =  40° 

H:  Huber  Solution  (Spherical  Indentation 
M:  Michell  Solution  (Point  Indentation) 


Figure  40.  Included  Angle  Effect  of  Conical  Indentation 

(I):  Non-dimensionalized  Stress,  a 
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Figure  41.  Included  Angle  Effect  of  Conical  Indentation 
(II):  Non-dimensionalized  Stress,  a . 
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Figure  42.  Included  Angle  Effect  of  Conical  Indentation 
(III) :  Non-dimensionalized  Stress,  o  . 
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Figure  43.  Included  Angle  Effect  of  Conical  Indentation 

(IV):  Non-dimensional  Stress,  a  . 
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2  2  2 
x  +  y  =  a  =  0.09 

By  comparing  these  figures  with  Figures  36  to  39,  it  is  obvious 
that  Sneddon's  stress  field  behaves  like  Hertzian  stress  field  under 
the  same  contact  radius  and  also  Sneddon's  stress  field  behaves  like 
Michell's  stress  field  under  the  sharp  included  angle  and  small 
contact  radius  as  be  expected. 
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APPENDIX  B 

THE  ANALYTICAL-EXPERIMENTAL  EQUATIONS  FOR  ELASTIC/PLASTIC  INDENTATION 


Various  formulas  to  define  the  indentation  hardness,  H,  in  terms 
of  Young's  modulus,  E,  Poisson’s  ratio,  v,  and  yield  strength,  Y, 
have  been  developed  by  Johnson  [11],  Perrott  [37,39],  and  Izumitani  [12]. 
These  formulas  were  derived  by  utilizing  Hill's  elastic-plastic 
solution  of  the  expansion  of  a  spherical  or  cylindrical  cavity  [7,8] 
along  with  experimental  observations. 

Johnson's  formula  for  a  conical  indentation  is: 


with 


Y  *  3  +  2  ln  (],■)  » 

|  tang  -  6(l-v)(^)3  -  4(l-2v)  , 


and  for  a  wedge  indentation  is: 


(B-l) 


(B-2) 


I  [\  +  m  (f)]  ,  (b-3) 

/3 

with 


H  tang  -  (5-4v)(|)2  -  3(1- 2v)  ,  (B-4) 

where  c  is  the  radius  of  the  plastic  zone  and  the  center  is  the 
undeformed  origin. 

Perrott 's  formula  for  conical  and  pyramidal  indentation  is: 


H 

Y 


0.494  +  0.577  in 


Etan 

Y(1-v2)J 


(B-5) 
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Izumitani's  formula  for  knoop  indentation  is: 


|  =*  j  +  2  In  (^)  , 


(£)3 _ S_ 

V  3(l-u)Y  * 


(B-6) 


(B-7) 


which  are  the  same  as  Hill's  equations  for  the  expansion  of  a  spherical 
cavity  in  an  infinite  medium. 

The  relationship  between  indentation  depth  vmax  and  residual  crater 
depth  vmaY  in  terms  of  material  properties  and  geometrical  properties 
for  the  conical  indentation  problem  is  expressed  by  Lawn  and  Howes 
[46]  as: 


-  “t3]l 


(B-8) 


where  H  is  the  indentation  hardness,  and  J£  is  the  elastic 
geometrical  parameter  which  is  expressed  as: 


E  btanB 


it  for  Sneddons  solution)!  , 


(B-9) 


and  Jg  is  actually  an  unknown  geometrical  parameter  of  plastic 


hysteresis;  however,  it  is  usually  defined  as  Y„  *  Y. 
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